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Abstract 

We investigate Bose-Einstein correlations (BEC) in W-pair production at ^J~s ~ 
189 GeV using the L3 detector at LEP. We observe BEC between particles from a 
single W decay in good agreement with those from a light-quark Z decay sample. 
We investigate their possible existence between particles coming from different Ws. 
No evidence for such inter-W BEC is found. 
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Introduction 



Bose-Einstein (BE) interference is observed as an enhanced production of identical bosons, 
e.g., charged pions, at small four-momentum difference in elementary particle and nuclear 
collisions [UI2], and, in particular, in hadronic Z decay [Hill]. Such an interference should 
also be present in hadronic W decay (intra-W BE interference). Furthermore, since in fully 
hadronic WW events (e + e~ — > W + W~ — > qqqq), the W decay products overlap in space-time, 
interference between identical bosons originating from different Ws can be expected jSHZ]- This 
inter-W BE interference may provide a laboratory to measure the space-time development of 
this overlap. Moreover, this effect, like colour reconnection [^1514121. can be a source of bias in 
the determination of the W mass in the four jet channel. Recent model predictions |HHZ|, as 
well as recent experimental results [TH] . are still contradictory. 

The main question we address in this paper is, therefore, whether inter-W BE interference 
exists. However, we also examine all BE interference, intra-W as well as inter-W, and make 
a comparison with that observed in hadronic Z decays, with and without the contribution of 
Z — > bb decays. 



Analysis Method 

Bose-Einstein interference manifests itself through correlations between identical bosons at small 
four-momentum difference. Correlations between two particles are described by the ratio of 
the two particle number density, p2(pi,P2), to the product of the two single particle number 
densities, Pi(pi)pi(p2)- Since we are only interested in Bose-Einstein correlations (BEC) here, 
the product of single particle densities is replaced by Po(Pi,P2), the two particle density that 
would occur in the absence of Bose-Einstein interference, resulting in the BE correlation function 

R2(Pl,P2) = —, 7 • (1) 

P0{Pl,P2) 

For identical bosons, R 2 — 1 is related to the space-time particle density through a Fourier 
transformation |2*|ll4j. 

Since we shall consider only pion pairs, the mass of the particles is fixed and the correlation 
function is defined in six-dimensional momentum space. Since Bose-Einstein correlations are 
largest at small four-momentum difference, Q = v/— (pi — P2) 2 , we parametrize R2 in terms of 
this single variable. While this is an oversimplification, as recent two- and three-dimensional 
analyses have shown [I], lack of statistics prevents such multi-dimensional analyses here. 

The following method [TH] is used to study inter-W BEC. If the two Ws decay independently, 
the two particle density in fully hadronic WW events, p 2 NW , is given by 

p] NW (PuP2) =pT{PiiP2)+pT{PiiP2)+pT{Pi)pT{P2)+pT{Pi)pT{P2) , (2) 

where the superscript, W + or W~, indicates the W which produced the particles. Assuming 
that the densities for W + and W~ are the same, eq. (J2J) becomes 

pl NW (PuP2) = 2pf( Pl , P 2) + 2pY( Pl ) pf(p 2 ) . (3) 

The terms p 2 vVW and p^ of eq. © are measured in the fully hadronic WW and the semi- 
hadronic events, respectively. To measure the product of the single particle densities, we use 
the two particle density p^bfQoi, P2) obtained by pairing particles originating from two different 
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semi-hadronic WW events (W + W~ — > £vqq), since by construction these pairs of particles are 
uncorrected. 

The hypothesis that the two Ws decay independently can be tested using eq. (J3J). In par- 
ticular, we write eq. (J3J) in terms of Q and use the test statistics 

Ap(Q) = pf w (Q) - 2pf(Q) - 2p^(Q) (4) 

and ww 

D(Q) = 2pW ( g) + 2pWW (Q) • ( 5 ) 

The advantage of this method is that it gives access to the inter-W correlations directly from 
the experimental data; there is no need for normalization by a Monte Carlo (MC) model. 

It is possible that the event mixing procedure introduces artificial distortions and that it 
does not fully account for some non-BE correlations or some detector effects. To diminish the 
effect of such inadequacies and to be able to compare more directly to other experiments, we 
also use the double ratio 

D\Q) = D D{Q \ 0) , (6) 

where Dmc, noBE is derived from a Monte Carlo sample with no BEC, or at least without inter-W 
BEC. 

In the absence of inter-W correlations, Ap = and D — D' — 1. To study BEC, we 
examine these relations for small values of Q, for like-sign particles. To judge the influence of 
other correlations on these quantities, we examine them also for unlike-sign particles and in 
Monte Carlo models. 



Data Selection 

The data used in this analysis were collected in 1998 by the L3 detector [TH], and correspond 
to an integrated luminosity of about 177 pb" 1 at a centre-of-mass energy of ^/s ~ 189 GeV. 

To obtain the two W + W~ event samples, one fully hadronic and the other semi-hadronic, 
we reconstruct the visible final state fermions, i.e., electrons, muons, r jets (corresponding to 
the visible r-decay products) and the hadronic jets corresponding to quarks, and apply the 
selection criteria described in Ref. with the additional requirement for the fully hadronic 
channel that the neural network output must be greater than 0.6. In total, 1032 semi-hadronic 
events and 1431 fully hadronic events are selected. 

The event generator koralw ^H] is used to simulate the signal processes. Within koralw 
BEC are simulated using the so-called BE 32 or BE algorithms [7]. For most comparisons 
we will use BE 32 , since it agrees better with the data. Where BE is used it will be explicitly 
stated. The BEC are implemented for all particles, which we refer to as BEA, or only for 
particles coming from the same W (intra-W BEC), which we refer to as BES. The background 
processes e + e~ — > Z/7 — > qq, e + e~ — > TIL and e + e~ — > Ze + e~ (the last relevant only to the 
qqeu and qqrz/ channels) are generated using pythia ^H] with BE . The generated events are 
passed through the L3 detector simulation program [20], reconstructed and subjected to the 
same WW selection criteria as the data. 

*)The BE32 algorithm used the parameter values PARj(92)=1.68 and parj(93)=0.38 GeV, whereas BE used 
PARj(92)=1.50 and parj(93)=0.33 GeV. Both have been tuned to the L3 Z decay data. 
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MC studies using the above generators show that the selection efficiency for fully hadronic 
events changes by less than 0.5% when BEC (intra-W, or both intra-W and inter-W) are in- 
cluded. The efficiencies for the channels qqiv, (£ = e,/z, r) and qqqq are found to be 82.9%, 
76.7%, 50.6% and 87.2%, respectively. The fractions of background for these channels are 4.1%, 
4.1%, 12.7% and 18.6%, respectively. 

The BEC study is based on charged particle information from the central tracker. Charged 
tracks are required to have at least 35 (of 62 possible) hits, and the number of wires from the 
first to the last hit is required to be at least 45. The distance of closest approach (projected 
onto the transverse plane) of a track to the nominal interaction vertex is required to be less 
than 7.5 mm. The transverse momentum of a track must be greater than 100 MeV. After the 
track selection, there are 287k pairs of like-sign particles in the fully hadronic channel and 55k 
pairs in the semi-hadronic channel. 

With this selection, reasonable agreement is obtained between the data and the MC sim- 
ulation for the distributions of Q and the difference in azimuthal, as well as polar, angle with 
respect to the beam, for pairs of like-sign tracks, in both the fully hadronic and the semi- 
hadronic channels. This is shown in Fig.Q where the raw data are compared to simulated 
koralw with BES and background events. Similar agreement is observed when BEA is used. 



Measurement of R2 

We first measure the BE correlation function, R2, for like-sign charged pion pairs using two 
choices of reference sample, i.e., the sample from which p is determined. The first choice uses 
a Monte Carlo model without BEC: 



Po(=t) =t) — P 2 (±) ±)MC,noBE 



(7) 



The second choice uses unlike-sign particle pairs from the experimental events. A major draw- 
back of this method is that the correlation function is affected by the presence of dynamical 
correlations, such as the decay of resonances. To compensate for this, the density for unlike-sign 
pairs is multiplied by the ratio of the densities for like- and unlike-sign pairs determined from 
a Monte Carlo model without BEC: 



Po(±, ±) = P 2 (+, -) 



P2(±,±) 
LP2( + ,-) 



(8) 



MC, noBE 



In both cases we need to correct the correlation function, R2, for detector resolution, ac- 
ceptance, efficiency and for particle misidentification. For this we use a multiplicative factor 
derived from Monte Carlo studies. Since we do not perform explicit hadron identification, this 
factor is given by the ratio of p 2 (±, ±) and p2(±, ±)/p 2 (+, — ), respectively, found from Monte 
Carlo, for pions at generator level to that found using all particles after full detector simulation, 
reconstruction and selection. Thus, using eqs. ((ZJ) and (jHJ) leads, respectively, to 
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Background is taken into account by replacing />2data in the above equations by 



1 ( dn dn hg \ 

P2data-b g - ^ ^ - -jq) , (11) 

where V is the purity of the selection, n is the number of pairs of tracks in the N ev data events, 
and nb g is the number of pairs of tracks corresponding to (1 — V)N ev background events. The 
background is estimated using Monte Carlo. In determining R 2 using eq. (|9|L we use koralw 
without BEC as the reference sample. For the detector correction, BES with the BE 32 algorithm 
is used for both eqs. © and (fTUj) . 

Fig.[2]shows the correlation function, eq. (jH}, for the fully hadronic and for the semi-hadronic 
WW events. We parametrize the Bose-Einstein enhancement at low Q values by 

R 2 {Q) =7(l + <$Q)(l + Aexp(-i? 2 Q 2 )) , (12) 

where 7 is an overall normalization factor, the term (1 + 5Q) takes into account possible long- 
range momentum correlations, A measures the strength of the BE correlations and R is related 
to the source size in space-time. The results of the fits of eq. (jl2J) are also shown in Fig-El 

The fit results for both choices of reference sample, eqs. (J7J) and (jHJ), are given in Tabled The 
statistical error includes bin-to-bin correlations. These are estimated from 100 sets of W + W~ 
BES events generated by pythia, each with the same statistics as the data. The variation of 
7, A, R and 5 from their average values was determined for the fully- and semi-hadronic WW 
events. The ratio of the Gaussian width to the average fit error was found to be between 1.01 
and 1.61. For each parameter, the corresponding ratio is used to scale the original statistical 
error. MC studies show that this ratio hardly depends on Q, which justifies this method to 
correct for bin-to-bin correlations. 

The systematic uncertainty is computed by varying the track and event selections. Both 
stronger and weaker cuts are applied to the tracks, slightly different event selections are made, 
and the background fractions are varied. The influence of the choice of the Monte Carlo used 
for the reference sample and for the correction factor is also taken into account. Part of the 
systematic uncertainty comes from the choice of the fit range. The large systematic uncertainty 
on A in the fully hadronic channel is mainly due to the difference of including or not including 
inter-W BEC in the MC for the correction factor. Using inter-W BEC in the correction factor 
increases the measured value of A. 

BEC are observed (A > 0) in both fully hadronic and semi-hadronic WW events. The values 
of A are higher for the semi-hadronic than for the fully hadronic channel, but the difference is 
only about two standard deviations (statistical error only) for each choice of reference sample. 
If true, this effect would indicate a suppression of inter-W BEC [T^], which we study in detail 
in the following section. 

Since, apart from the quark flavour, hadronic W and Z decays are expected to be similar, we 
also analyze a high statistics hadronic Z decay sample, collected by the L3 detector in 1994 at 
y/s ~ 91.2 GeV. Since b quarks are greatly suppressed in W decays, a b-tagging procedure [21] 
is used to reduce the bb fraction in Z decays, from 22% to 3%. The BE correlation function, 
eq. (JHJ), of the resulting 180k Z events is plotted in Fig.Eb as a full histogram. As expected, good 
agreement is observed between this histogram and the correlation function of the semi-hadronic 
WW events. When b quark decays of the Z are not removed from the sample, a depletion of 
the correlation function at small Q is observed and a clear discrepancy exists with the W data, 
as the dashed histogram in Fig.Eb shows. 
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Measurement of Inter-W Bose-Einstein Correlations 



The Event Mixing Procedure 

To compute the test statistics, eqs. (jlj and ©, we need to construct the two particle density 
Pmh? ■ This is done by combining pairs of semi-hadronic events having oppositely charged 
hadronically decaying W's. Particles identified as decay products of the leptonically decaying 
W's are discarded. Then the particles from one of the events are rotated so that the W's are 
approximately back-to-back. Since real fully hadronic WW events have a small longitudinal 
energy imbalance that we ascribe to initial state radiation and since experimental resolution 
leads to both transverse and longitudinal energy imbalance, we do not force the Ws to be 
exactly back-to-back. We introduce an extra momentum, p ext r&, Gaussian distributed in all 
three components and impose p ex tra + W\ = —W2, where Wi^ are the momenta of the two 
Ws. For the longitudinal component the Gaussian has mean and standard deviation 7.9 
GeV, while for the transverse components the mean is randomly chosen as ±0.5 GeV and 
the standard deviation is 1.4 GeV. These values were chosen to obtain reasonable agreement 
between the energy imbalance distributions of fully hadronic and mixed events. 

In addition, we impose the following cuts which are related to the pre-selection of fully 
hadronic WW events ^7j. We demand that the sphericity be larger than 0.045, that the total 
visible energy be larger than 0.7 y/s, that the number of particles identified with the calorimeter 
(the cluster multiplicity) be larger than 30, that the ratio of the total longitudinal energy 
imbalance to the visible energy be smaller than 0.25, and that the y cut value at which the event 
changes from a 3- to a 4-jet topology, 2/34, be larger than 0.001. After forcing the event into 2 
jets with the Durham clustering algorithm [22], the average of the jet masses is required to be 
larger than 30 GeV. After forcing the event into 4 jets with the Durham clustering algorithm, 
we assign two pairs of jets to each of the two Ws by first rejecting the combination with the 
smallest dijet mass and then accepting the combination with the smallest difference between 
the two dijet masses (best pairing). We then demand that the difference between the two W 
masses be less than 70 GeV, that the smallest angle between any two jets be larger than 0.28 
radians, and that the average of the two smallest angles between two jets from different Ws 
be larger than 0.6 radians. These cuts reject only approximately 1% of the events and do not 
change the Q-distribution. The final selection of fully hadronic events uses a cut at 0.6 on the 
output of a neural network [T^j. This cut is also applied to the mixed events, rejecting 7%. 

We have checked the mixing procedure by comparing the distributions and quantities of a 
large number of variables between mixed events and fully hadronic WW events, including event 
shape variables, track and cluster multiplicities, and variables related to the W such as mass, 
energy and orientation. In general, good agreement is found. Typical examples are shown in 
Fig.El 

Results 

Fig.0] shows the distributions of the three terms in the right-hand side of eq. (j3J) for the data, 
not corrected for detector effects, but after background subtraction. At low values of Q we 
observe more pairs of unlike-sign particles than pairs of like-sign particles, both in the two- 
particle densities for fully hadronic (Fig.El) and semi-hadronic (Fig.Hb) events. Furthermore, 
we observe that p^(±, ±) and p^^(+, — ) coincide (Fig.Efc). 

From these distributions, we compute Ap for like-sign and unlike-sign particle pairs, eq. (J4*j). 
The resulting raw data distributions are shown in Fig. 03 Also shown are the predictions of 
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koralw after full detector simulation, reconstruction and selection. Both the BEA and BES 
scenarios are shown. 

The BEA scenario using BE32 shows an enhancement in the Ap distribution for like-sign 
pairs (Fig.Et)j but also a small enhancement for unlike-sign pairs (Fig.Eb)- The effect for 
unlike-sign pairs is larger if BEo is used. These implementations of BEC clearly affect both 
the like- and unlike-sign particle spectra. From Fig.^ it is clear that only the BES scenario 
describes the Ap(±, ±) distribution, while the BEA scenario is disfavoured. In particular, 
considering Q-values up to 0.6 GeV, the confidence level (CL) for the BES and BEA scenarios 
are, respectively, 84% and 0.8%. The calculations of the confidence levels are based on statistical 
errors including bin-to-bin correlations. 

In Fig.|Hlwe show the distributions of D and D' for like-sign and unlike-sign particle pairs, 
eqs. (JHJ) and (J0J), for the raw data. For the double ratio D' we use the BES scenario of koralw 
as the reference sample. Also shown in the figure are the predictions of koralw for the 
scenarios BEA and BES. Again, it is clear that the BES scenario of koralw describes the data, 
CL=87% for both D(±,±) and D'(±,±), while the BEA scenario is disfavoured, CL=0.5% 
for both D(±, ±) and £)'(±, ±). When BEo is used instead of BE32, the BEA scenario is 
even more strongly disfavoured: CL=0.08% for both D(±, ±) and -D'(±, ±). Note that the D' 
distributions are by definition equal to unity (apart from statistical fluctuations) when koralw 
without inter-W BEC is used. Note also that D is already close to unity for BES, so that the 
difference between D and D' is small, which supports the validity of the mixing procedure. 

To estimate the strength of inter-W BEC, the -D'(±, ±) distribution is fitted (from to 1.4 
GeV) by the following function 

D'(Q) = (1 + eQ) (l±Aexp(-A; 2 Q 2 )) , (13) 

where e, A and k are the fit parameters. The result of the fit for the strength of inter-W BEC 
is 

A = 0.001 ±0.026 ±0.015 , 

where the first error is statistical and the second systematic. The statistical error has been 
multiplied by 1.49 to account for bin-to-bin correlations, in the same way as described in the 
previous section. This value of A is consistent with zero, i.e., with no inter-W BEC. A similar 
fit was performed for the koralw BEA distribution, resulting in A = 0.127 ±0.007 (statistical 
error only). The data disagree with this value by more than 4 standard deviations. 

The systematic uncertainty on A is the sum in quadrature of the contributions listed in 
Table 121 The amount of background was varied by ±10%. The choice of Monte Carlo was 
varied using pythia and koralw, both with no BEC at all as well as with only intra-W 
BEC. Also the effect of various models of colour reconnectiorjl3 (CR) was included. A change 
in the fit range (±0.4 GeV), a change in the bin size (from 40 to 80MeV) and a change in 
the parametrization (removing the factor (1 ± eQ) from the fit) also give contributions to the 
systematic uncertainty. Furthermore, the track and event selections were varied. 

In the mixing procedure we allow a semi-hadronic WW event to be combined with all 
possible other semi-hadronic WW events. To be sure that this does not introduce a bias, 
the analysis was repeated for a mixed sample where every semi-hadronic event was used only 
once. The influence of the mixing procedure was also studied by not only combining oppositely 
charged Ws, but also like-sign Ws. The influence of the extra momentum p ex tra, used in the 
event mixing, is also included as a systematic effect. The RMS of the systematic uncertainties 

t)The so-called SKI, SKII, SKII' and GH (TH] models, as implemented in PYTHIA, were used. 
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due to these three changes in the mixing procedure is the systematic uncertainty listed in 
Table El The influence of the cut on the neural network output for the mixed events was 
investigated by removing the cut. 

Furthermore, the effect of uncertainties in the energy calibration of the calorimeters was 
studied. Finally, we studied the influence of the qqr^ channel. Since this channel is the most 
difficult to identify, and therefore has relatively high background and low efficiency, we repeated 
the analysis without it. 

To make the analysis possibly more sensitive to inter-W BEC, we repeated the analysis twice 
using different selections to increase the overlap of the W + and W~ decay products. Since BEC 
occur mainly among soft particles and the overlap is expected to be larger for these particles 
than for high-momentum ones, we first repeated the analysis using only tracks with momenta 
smaller than 1.5 GeV. Another way to increase the overlap is to require that jets from different 
W's be close together. We therefore repeated the analysis requiring that the average of the 
smallest two of the four angles between jets from different W's be less than 75°. This results 
in a reduction of approximately 60% in the number of fully hadronic WW events. Fig. [7| shows 
the distributions of D'(±, ±) for these two analyses. It is again clear that the BEA scenario is 
disfavoured, particularly for the low-momentum selection, while BES describes the data well. 
For the low-momentum sample we find CL=1.6% for BEA and CL=96% for BES, and for the 
sample with the angular cut we find CL= 10% for BEA and CL=93% for BES. Moreover, we 
find A = 0.026 ± 0.034 for the low-momentum sample and A = —0.019 ± 0.029 for the sample 
with the angular cut. Both values are consistent with zero. The errors here are statistical only, 
including bin-to-bin correlations. 

Conclusion 

Intra-W Bose-Einstein correlations have been found to be similar to those observed in Z decay 
to light quarks. An excess at small values of Q in the distributions of Ap(±, ±), D(±, ±) and 
-D'(±, ±) is expected from inter-W BEC, but none is seen. These distributions agree well with 
koralw using BE 32 when inter-W BEC are not included, but not when they are. We thus find 
no evidence for BEC between identical pions originating from different Ws and disfavour their 
implementation using the BE 32 and BE algorithms. 
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reference 


MC, no BEC 


(±, -) pairs 


channel 


fully-hadronic 


semi-hadronic 


fully-hadronic 


semi-hadronic 


7 
A 

R (fm) 
5 

xVriclf 


0.91 ±0.02 ±0.02 
0.55 ±0.04 ±0.08 
0.56 ±0.04 ±0.06 
0.06 ±0.02 ±0.06 
32/31 


0.94 ±0.01 ±0.02 
0.70 ±0.06 ±0.05 
0.64 ±0.05 ±0.06 
-0.01 ±0.01 ±0.06 
41/31 


0.93 ±0.01 ±0.02 
0.48 ±0.05 ±0.08 
0.71 ±0.04 ±0.05 
0.07 ±0.01 ±0.05 
29/31 


0.89 ±0.01 ±0.03 
0.64 ±0.07 ±0.06 
0.75 ±0.05 ±0.06 
0.07 ±0.02 ±0.05 
35/31 



Table 1: Values of the fit parameters 7, A, R and 5, eq. (JUJ), for the fully hadronic and 
the semi-hadronic WW events. Two different reference samples are used: koralw 
without BEC, eq. (jZJ), and unlike-sign particle pairs, eq. (jHJ). The first error is 
statistical, the second systematic. 



source 


contribution 


background fraction 


0.0021 


other Monte Carlo (pythia, BES or no BE) 


0.0060 


allowing CR in the reference sample 


0.0024 


fit range 


0.0012 


rebinning (40 -> 80 MeV) 


0.0024 


removing (1 ± eQ) from the fit 


0.0011 


track selection 


0.0073 


event selection 


0.0046 


mixing 


0.0040 


neural net output cut 


0.0039 


energy calibration 


0.0014 


influence of r channel 


0.0076 


total systematic uncertainty 


0.015 



Table 2: Contributions to the systematic uncertainty of the A parameter. Explana- 
tion of the sources are in the text. 
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Figure 1: Distributions of (a, d) 5<p, the difference in azimuthal angle of pairs of 
tracks, (b, e) 89, the difference in polar angle of pairs of tracks, and (c, f) Q, the 
four-momentum difference, for the fully hadronic WW events (a-c) and for the 
semi-hadronic WW events (d-f). Only like-sign pairs of tracks are considered. The 
points are the uncorrected data, the open histograms are the expectation of KOR- 
alw with intra-W BEC plus background. The shaded histogram is the background 
expectation. 
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Figure 2: The Bose-Einstein correlation function R 2 , eq. ({£)]). for (a) the fully- 
hadronic WW events, and (b) the semi-hadronic WW events. In (b) the full his- 
togram is for the light-quark Z decay sample and the dashed histogram is for a 
sample containing all hadronic Z decays. Also shown are the fits of eq. (jl2j) to the 
WW data. 
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Figure 3: Comparison of uncorrected distributions for fully hadronic events after 
background subtraction (points) and mixed events (histograms): (a) — log 1/34; (b) 
the two smallest angles between jets of different W's, after jet finding and best 
pairing; (c) the cosine of the angle ip between the decay planes of the two W's, after 
jet finding and best pairing; and (d) the event thrust. 
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Figure 4: Distributions for uncorrected data of (a) p^ w , (b) and (c) for 
pairs of like-sign charged particles and pairs of unlike-sign charged particles. 
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Figure 5: Distributions for uncorrected data of (a) Ap(±, ±) and (b) Ap(+, — ). 
Also shown are the Monte Carlo predictions of koralw (at the detector level) with 
BEA (inter-W) and BES (no inter-W). 
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Figure 6: Distributions for uncorrected data of (a) D(±,±), (b) D(+,—), (c) 
D'(±, ±) and (d) D'(+,—). Also shown are the predictions of koralw (at the 
detector level) with BEA (inter-W) and BES (no inter-W). 
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Figure 7: Distributions for uncorrected data of D'(±, ±) where (a) only low momen- 
tum tracks are used and (b) a cut is made on the average angle of the two smallest 
angles between jets of different W's. Also shown are the predictions of koralw (at 
the detector level) with BEA (inter-W) and BES (no inter-W). 
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